Excited States Calculated by Means of the 
Linear Muffin-Tin Orbital Method 

M. Alouani 1 and J. M. Wills 2 

1 IPCMS, Universite Louis Pasteur, 23 Rue du Loess, 67037 Strasbourg, France 

2 Los Alamos National Laboratory, Los Alamos, NM 87545, USA 



Abstract. The most popular electronic structure method, the linear muffin-tin 
orbital method (LMTO), in its full-potential (FP) and relativistic forms has been 
extended to calculate the spectroscopic properties of materials form first principles, 
i.e, optical spectra, x-ray magnetic circular dichroism (XMCD) and magneto-optical 
kerr effect (MOKE) . The paper describes an overview of the FP-LMTO basis set and 
the calculation of the momentum matrix elements. Some applications concerning 
the computation of optical properties of semiconductors and XMCD spectra of 
transition metal alloys are reviewed. 



1 Introduction 

The density functional theory (DFT) of Hohenberg, Kohn, and Sham is the 
method of choice for describing the ground-state properties of materials H. 
However, in the initial derivation of the DFT, the eigenvalues are Lagrange 
multipliers introduced to orthogonalize the eigenvectors, which in their turn 
are used to compute the total energy and the charge density. In this formu- 
lation the eigenvalues have therefore no physical meaning and should not be 
considered as excited states. Nevertheless, the DFT in the local density ap- 
proximation (LDA) or in its spin resolved local density formulation (LSD A), 
has been used successfully to compute the excited states, namely, optical and 
magneto-optical properties, x-ray absorption and magnetic dichroism spec- 
tra. 

The LDA or LSDA were indeed intended to compute the ground-state 
properties of materials, and their use during the last two decades has pro- 
duced an excellent track record in the computation of these properties for a 
wide variety of materials, ranging from simple metals to complex semicon- 
ductor superlattices. However, it is now believed that the DFT can do more 
than computing the ground state properties. This is because the Kohn-Sham 
equations could be viewed as deriving form a simplified quasi-particle (QP) 
theory where the self-energy is local and time averaged, i.e., S(r,r',t) « 
V xc (r)5(r — r')5(t), here V xc (r) is the local exchange and correlation poten- 
tial as, for example, parameterized by Von Barth and Hedin Viewed in 
this way, the KS eigenvalues are then approximate QP energies and could be 
compared to experimental data. This argument is supported by quasiparticle 
calculations within the so called GW approximation of Hedin [|| showing 
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that the valence QP energies of semiconductors are in good agreement with 
these obtained using LDA, and the conduction QP energies differ by approx- 
imately a rigid energy shift In the literature this shift is often called 
"scissors-operator" shift Q. 

In the last few years spectroscopy is becoming the standard tool for mea- 
suring excited states of materials. Its owes its impressive advances to the 
availability of synchrotron tunable highly polarized radiation. In particular, 
the measurement of optical, magneto-optical properties as well as magnetic x- 
ray dichroism are now becoming routine tasks for probing the structural and 
magnetic properties of materials. Considerable attention has been focused 
on transition-metal surfaces and and thin films due to their novel physical 
properties different from that of bulk materials and due to potential indus- 
trial applications such as magneto-optical recording, sensors, or technology 
based on giant magneto-resistance. In this respect, theory is falling far be- 
hind experiment and it is becoming hard to give a basic interpretation of 
experimental data. 

This paper, which is far from being a review paper about calculated ex- 
cited states, tries to bridge the gap between experiment and theory by describ- 
ing a rather quantitative method for computing excited states of materials. 
This method uses the local density approximation and the linear muffin-tin 
orbital (LMTO) method. In the first part of this paper we introduce the 
density functional theory and the local density approximation and justify 
the use of LDA eigenvalues as approximate excited states and relate them 
to quasiparticle energies. In the second part we describe the construction of 
the LMTO basis set within an all-electron full-potential approach which 
will be used to determine the momentum matrix elements. We devote the 
third part to the determination of the momentum matrix elements. In the 
first part of the application section we present some examples of computation 
of semiconductors optical spectra [p|jic| , and leave out the optical properties 
of metals and magneto-optical properties of materials and refer the reader 
to Ref. Jll|,||Jl|@|^|60,|l|,|l9|. In the second part of the applications we 
show some examples of x-ray magnetic dichroism calculations |2(],|2l],^2| . 



2 Density Functional Theory 



The density functional method of Hohenberg and Kohn |lj which states that 
the ground state total energy of a system of N interacting electrons in an 
external potential V^ xt is a functional of the electron density p(r) does not 
provide an analytical form of the functional [Q . This method remains numer- 
ically intractable without the Kohn and Sham introduction of the so called 
local density approximation [|l| in which the exchange and correlation func- 
tional E xc {n} appearing in the total energy: 
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E{n} = T{n} + ^ J d\ J dV^^ + E xc {n} 

d 3 r V cxt (r)n(r) + E(V cxt ) (1) 



is given by E xc {n} — J d 3 r e xc (n(r))n(r) where e xc is the exchange- 
correlation energy of a uniform electron gas of density n. Thus, Kohn and 
Sham constructed a set of self-consistent single-particle equations: 

-V 2 + ^(£-T)W) =e^(r) (2) 



5n 

where the density n(r) is given by: 

n{v) =Y,0{ti < E F )U*)4{?) (3) 



and 



Rt 



E(V cxt ) = e 2 £ ]T (l - R ') S ^ O) | T + R r /; , _ R1 © 

Instead of the true kinetic energy of the electron gas, Kohn and Sham 
used the homogeneous electron kinetic energy: 

f = Y,0(e*<E F ) fd\4(-V 2 )^ (6) 

This use of homogeneous-electron kinetic energy in the Kohn-Sham equa- 
tions redefined the exchange-correlation function to be: 

E xc {n} = (E xc {n} + T{n} - f) = J d 3 r e xc {n{v))n{v) (7) 

It is then crucial to use a good basis-set for the description of the elec- 
tronic structure of realistic systems. The augmented plane wave Q (APW), 
and the Korringa-Kohn-Rostoker |24j] (KKR) methods can be used, in prin- 
ciple, to solve exactly the Kohn-Sham equations, however these methods are 
numerically involved and their linearization, introduced by Andersen is much 
preferable. Andersen linearization, has not only made the techniques for solv- 
ing the band-structure problem transparent by reducing it essentially to the 
diagonalization of one-electron Hamiltonian, and cuts the cost of computa- 
tion by at least one order of magnitude. The linearized versions of these two 
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powerful methods are the linear augmented plane wave (LAPW) and linear 
muffin-tin orbital (LMTO) methods, respectively |||. 

In this paper, we will only use the LMTO method to study excited states 
of solids. The reason for this choice is that the LMTO method is the mostly 
used method in computational electronic structure. This is due primarily 
to the use of atomic-sphere approximation (ASA) which made the LMTO 
method run fast even on today's cheap personal computers. Due to this re- 
duced computational cost, the LMTO ASA method became the method of 
choice of researchers without access to supercomputers. 



3 Quasiparticle Theory and Local-Density 
Approximation Link 

The quasiparticle (QP) electronic structure of an interacting many-body sys- 
tem is described by the single-particle eigenstates resulting from the interac- 
tion of this single particle with the many-body electron gas of the system. The 
single-particle eigenstate energies are the results of solving a Schroedinger like 
equation containing the non-local and energy-dependent self-energy instead 
of the exchange-correlation potential appearing in Kohn-Sham like equations: 

(T + V H + V c ^{v) + J d 3 r'S(r, r', E)V{r') = m{r). 

Thus the self-energy 27 contains all many-body effects. Almost all ab- 
initio QP studies were performed within the so-called GW approximation, 
where the self-energy 27 is calculated within Hcdin's GW approximation. This 
method consists of approximating the self-energy as the convolution of the 
LDA self-consistent Green function G and the screened coulomb interaction 
W within the random-phase approximation. The QP eigenvalues are often 
obtained using first-order perturbation theory starting from LDA eigenvalues 
and eigenvectors [p|]25||. Although there are early calculations starting from 
Hartree-Fock [|(| or tight-binding p7f methods. Nevertheless, the best results 
are based on a LDA starting point [[l[ |25| , |28| , p9 30 . 



Thus the GW predicted optical excitations energies of semiconductors are 
within 0.1 eV form the experimental results and the surprizing fact is that 
the QP wave functions are almost identical to these produced within the 
LDA H (the wave function overlap is more than 99%). For a general review 
of GW calculations see the review by Araysetianwan and Gunnarsson [^9| or 
by Aulbur, Jonsson and Wilkins [£30| . 

It is clear that the quasiparticle Schroedinger equation resembles to the 
Kohn-Sham equation. Both equations describe a fictitious electron moving 
in a effective potential. The difference is that the self-energy is nonlocal and 
energy dependent whereas the LDA potential is local and averaged over time. 
This resemblance can be further pushed by noticing that the DFT can be used 
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to obtain excitation energies. For example, the ionization energy, /, and the 
electron affinity, A, are difference between ground state energies: 

I = E(N-1)-E(N), and A = E(N) - E(N + 1) 

Were N is the number of electrons of the system. And since the DFT 
gives the correct ground state energies it should produce, in principle, the 
correct ionization and electron affinity energies. For metals, the addition or 
removal of an electron from the system costs the same energy, and hence the 
ionization energy is equal to the electron affinity. For insulators, the energy 
gap makes all the difference and hence breaks this symmetry. Thus the energy 
band gap is given by: 

Eg = I — A = E(N + 1) + E(N - 1) - 2E(N) 

In practice, however, the calculation is often obtained within the LDA 
and the energy band gap is calculated as the difference between the lowest 
conduction band and the highest valence band. It was shown by Sham and 
Schluter |3l[] and Perdew and Levy |}2| that the calculated energy gap differ 
from the true band gap by an amount A even when the DFT is used without 
the LDA. The A value could range from 50% in the case of silicon to 100% in 
the case of germanium. For most of the semiconductors, the GW calculations 
show that the LDA eigenvalues differ form the GW quasiparticle energy by 
a constant A which is almost independent of the k-point. This finding is 
important and shows that the LDA eigenvalues have some meaning and could 
be used to calculate excited states. So as stated in the introduction, the 
Kohn-Sham equations could be viewed as deriving form a simplified quasi- 
particle (QP) theory where the self-energy is made local and time averaged, 
i.e., £(r,r',i) ~ V xc (r)S(r — r')6(t). This approximation is certainly good for 
metals where we have a good data base for excited state calculated within the 
LDA 0|||jl|Q||,|6|0,|l|,|§ and where the agreement with experiment is 
good. For semiconductors, this approximation is not bad either, provided we 
know the value of the discontinuity of the exchange and correlation. Usually, 
this value is provided by GW calculations or by experiment. 

4 The Full-Potential LMTO Basis Set 

In this section we describe the LMTO basis-set used to calculate the excited 
states of solids. We discuss the basis used for an all electron calculation where 
the potential is not supposed to be spherically symmetric nor of muffin-tin 
type. The use of a general potential makes the study of open structures possi- 
ble without having to resort to the so-called "empty-sphere" approximation. 
To define the basis-set, we divide the space into non overlapping spheres 
called "muffin-tin" spheres and a region between these spheres which we call 
interstitial region. Inside the muffin-tin spheres the Schroedinger equation is 
solved at a fixed energy for each angular momentum t and variational pa- 
rameter k (which is defined later). The linearization amounts to the use of a 



6 



M. Alouani and J. M. Wills 



linear combination of the solution <f)e(e, r) of the Schroedinger equation for a 
fixed energy and its energy derivative <j)e(e,r) inside the muffin-tin spheres. 
These linear combinations matche continuously and differentiably to an en- 
velop function (spherical function) in the interstitial region. The Bloch wave 
function in the interstitial region is given by a linear combination of these 
Hankel functions centered at each site: 



&(k, r) = e ik - R ^L, P-ri-R) (8) 

R 

where i stands for the number of the basis function quantum numbers (these 
numbers are {r, L, k, {ett}}), where r is the site number, L — (£,m) groups 
the two angular quantum numbers, and en is the linearization energy for 
a particular atom type t and angular momentum number I. The envelop 
functions are defined as /Cl(k, r) = Kt (k, r)IVi(r). 



y tm (v) = m m {v) (9) 

K t {n,r) ^ -K l+1 \ m{ r\ -r , ?A >0 n < i h (10) 
y ' ; [ne(Kr) - ije(Kr), if k j < 0, (« = i\k\) y ' 

Jt(K,r) = K-'Mitr) (11) 

Here ni is the Neumann function and ji Bessel function for the angular 
momentum £, and Yg m are the spherical harmonics. 

To get the differentiability of the wave-function at the boundary of the 
muffin-tin spheres, we write the envelope function inside the muffin-tin spheres. 
The envelope function for a muffin-tin sphere r' is given by: 



V e lkR /C L ( K ,r-r-R) (12) 
r 

= ]T y L >{T T ,) (jCt>(K, r T )S(r, t')S(L, L') (13) 
V 

+ J#(k, r T )B L f >L (T'-T, k, k) 

To produce smooth basis functions we require that the basis function is 
differentiable at the boundary of each muffin-tin sphere, i.e., that a linear 
combination of <f> and <f> matches continuously and differentiably K, and J 
at the boundary of the parent sphere and other spheres, respectively. Using 
these matching conditions at the muffin-tin spheres, the Bloch wave function 
inside a muffin-tin sphere r of the unit cell at the origin is given by : 

V>i(k, r) = J~] yL{V T r T )Ui>{eiu-, r T )fi(lt, e a iKi)B L ,L i ( T ~ T i' K ^ k )( 14 ) 
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where 

^(e,r) = (^(e,r),0/(e,r)) (15) 

To add the spin dependence to the basis-set, the Bloch wave function is 
multiplied by the eigenvector of the Pauli spin operator rj±i: 

Vv(k,r) = ^(k,r)r? CT 

Such that: 

n • 077 ±a = (±l)?7 ± i 
where 77 is the quantization axis chosen in advance. 



5 Dielectric function 



5.1 Dynamical Dielectric Function 

Here we give a review of the determination of the dielectric response of a 
semiconductor due to the application of an electric field. We expend the 
description of our published work []Tc| by giving more details concerning the 
calculation of the momentum matrix elements. 

An electromagnetic field of frequency u>, and a wave vector q + G inter- 
acting with atoms in a crystal produces a response of frequency uj and a wave 
vector q + G' (G and G' being reciprocal lattice vectors). The microscopic 
field of wave vector q + G' is produced by the umklapp processes as a result 
of the applied field E (c[ + G, uj) 

E {q + G,uj) =^e G ,G'(q,w)^(q + G',c J ) (18) 

G' 

where E(q + G, ui) is the total field producing the non-diagonal elements 
in the microscopic dielectric function £g,G'(ci, w). The microscopic dielectric 
function in the random phase approximation is given by |33| : 



e G ,G'(q,w) = Sg,G> ~ 

x E 



k.n.n' 



8TTe 2 
^|q+G||q + G'| 

/n',k+q /n.k 
-En',k+q — ^n,k — Tu*) + %8 



(19) 



(n', k + q|e i ( q+G ) r |7i, k>(», k|e- j(q+G ' )r |?i', k + q) 
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Here n and n' are the band indexes, f n ^ is the zero temperature Fermi 
distribution, and fl is the cell volume. The energies E n ^ and the the crystal 
wave function |n,k) are produced for each band index n and for each wave 
vector k in the Brillouin zone. 

The macroscopic dielectric function in the infinite wave length limit is 
given by the inversion of the microscopic dielectric function: 



e(uj) = lim 



[ £ G !G'(q^)]o,o 
eo,oM - lim ^ eo,c(q, w)T~ 1 G ,(q, w)e G ',o(q, £ 

q ^ G,G'^0 



(20) 



Where Tqq, is the inverse matrix of Tg,g' containing the elements €g,g' 
with G and G' ^ 0. The first term of this equation is the interband contri- 
bution to the macroscopic dielectric function and the second term represent 
the local-field correction to e. The most recent ab-initio pseudopotential cal- 
culation found that the local-field effect reduces the static dielectric function 
by at most 5% [p|. Previous calculations with the same method have also 
found a decrease of Eoo by about the same percentage 0,[34| . For insulators 
the dipole approximation of the imaginary part of the first term of equation 
© is given by ||: 

e 2 f 

= ^i^X! / dk K n ' k l V l n '' k )| 2 /«,k(l-/n'.k)^(ek,n',«-^) ,(21) 

Here v is the velocity operator, and in the LDA v = p/m (p being the 
momentum operator), and where ek,n,n' = E n '^ — E n ^. The matrix elements 
(nk|p|n'k) are calculated for each projection pj = jdj, j = x or y and 
z, with the wave function |nk > expressed in terms of the full-potential 
LMTO crystal wave function described by equations (|l4| ) and (||). The k- 
space integration is performed using the tetrahedron method |36| with a large 
number of irreducible k points the Brillouin zone. The irreducible k-points 
are obtained from a shifted k-space grid from the high symmetry planes and 
-T point by a half step in each of the k x , k y , and k z directions. This scheme 
produces highly accurate integration in the Brillouin zone by avoiding high 
symmetry points. 



5.2 Momentum Matrix Elements 

To calculate these matrix elements we first defined a tensor operator of order 
one out of the momentum operator Vo = V 2 = J-j and V±i = T-^j( ± 

The muffin-tin part of the momentum matrix elements is calculated 
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using the commutator [V 2 ,^^] = 2V^ so that: 

drcj) Tll (r) Y t , ml (r^~T)X7 li (j) Tl (r)Yi> m (Y^T) = -fG^ l/ , m , 

where £ , m , are the usual Gaunt coefficients, and S r is the radius of 
the muffin-tin sphere of atom r. In the interstitial region the plane-wave 
representation of the wave function (see equation |J) makes the calculation 
straightforward, but a special care has to be taken for the removal of the 
extra contribution in the muffin-tin spheres. However, we find it much eas- 
ier and faster to transform the interstitial matrix elements as an integral 
over the surface of the muffin-tin spheres using the commutation relation of 
the momentum operator and the Hamiltonian in the interstitial region. The 
calculation of the interstitial momentum matrix elements is then similar to 
the calculation of the interstitial overlap matrix elements. The k — case 
has been already derived by Chen using the Korringa, Kohn and Rostokcr 
Greens- function method [B7| . We have tested that both the plane- wave sum- 
mation and the surface integration provide the same results. 

-V 2 = K 2 ptp 

A Hankel function can be integrated over a volume by knowing its integral 
over the bounding surface: 



= («? - 4) J cPii/faih (23) 

The surface of the interstitial consists of the exterior of the muffin-tin 
spheres and the unit cell boundary. 

Over the surface of the muffin tins: the surface area is S 2 df2 and the 
normal to the sphere points inward 

(n\ - 4) J d 3 r ^*p^ 2 = (24) 

At a muffin-tin sphere boundary S T the Bloch wave function is given by: 
^(k,r)| =^e ik - R /C Ll (^,r-T 4 -R) (25) 
= ^ yt™(r)~Ki(Ki, S)B emA „ H (r-Ti, Ki, k) 
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where B^^.^.k) = f^)^)«M)' 

Let W denotes the Wronskian W(f, g) = fg' — f'g 
We define then 



and K 



S 2 Wo = S 2 W(Kl(K)tK e ( K )) = 



1 

-1 



and 



w T eKi, k 2 = S T 



W(KJ( K )),K £ (k)-W 



K l K 2 



P^\r = E^ E K ^-lm-M ( o' 1 ) m_ '" ; ^' m; lj ^) 



K£+i m - M ( i ^2 ) G(t+1, m-fi; t, m; 1, ^) 



k; 



then 



EE^E 

r n im 
K? 



B^_i m _^ /Tn/ (r-T/, Kf,\s.)w T e-!(Kf, «i) 



x ( o iJ^(^- 1 . TO -w m ; 1» A*) 

-B^ + i m _^ /m/ (r-r/, k/, k)to T ^+i(K/, k,) 

X (o K?) m_ ^' ^' m ' 1; M ) 

B< mi < imi (T-Ti,Ki,k)^ +Z\(/,i,Kj) 

where 



+ B <,+im,- )1 i / m / (n-T/,K/,k)e(l 1 +l, m r /i; ^, m»; 1, /i)k 2 
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+ X> M (r / ) (31) 

+Bi f+ x nH + M ,i f m f (Tf-Ti,Ki,'k)G(tf, rrn\ £ f +l, mf+fi; 1, (j,) 
-B^_i m/ _ AlAmi (T/-r i ,/€ i ,k)^(^/ ! to/; £f-l, rrif+n; 1, ^)k 2 

5.3 Velocity Operator and Sum Rules 

Equation (|2l]) can not be used directly to determine the optical properties 
of semiconductors, when the GW approximation or the scissors operator is 
used to determine the electronic structure. The velocity operator should be 
obtained from the effective momentum operator p eff which is calculated using 
the self-energy operator, S(r, p), of the system Q: 

v = p efi 7TO = p/TO + <9Z , (r,p)/dp (32) 

GW calculations show that the quasiparticle wave function is almost equals 
to the LDA wave function Based on this assumption, it can be easily 

shown |38| that in the case of the scissors operator, where all the empty 
states are shifted rigidly by a constant energy A, the imaginary part of the 
dielectric function is a simple energy shift of the LDA dielectric function 
towards the high energies by an amount A, i.e., e^ P (cj) = e\ T:iA (u — A/fi). 
The real part of the dielectric function is then obtained from the shifted 62 
using Kramers-Kronig relations. The expression of e^ p is given by: 

qp 2e 2 ^ f |(n,k|pK,k)| 2 

& =1 + ^vEj */n.k(l - M ( + A)&1 , (33) 

n,n' ■ K -'' t ) ft 

e^ p is very similar to e^P A except that one of the interband gap Ck,n',n is 
substituted by the QP interband gap ek, n ',n + A. 

To test for the accuracy of the calculation within the LDA the f-sum rule: 

^EE - /n^) Kn ' k|pKk)|2 = 1, (34) 

k n,n 

where n v is the number of valence bands, should be always checked to ensure 
the accuracy of the calculations. 

It is easily seen that the dielectric function e% calculated using the 
scissors-operator shift does not satisfy the sum rule (up is the free-electron 
plasmon frequency): 

ue2(u)du = -^°p (35) 

because (i) ey° A satisfies this rule, and (ii) eJf is obtained by a simple shift 
of €2 DA by the scissors-operator A towards higher energies. The non simul- 
taneous satisfaction of both the f-sum rule and the integral sum rule within 
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the scissors approximation shows the limitation of this approximation. While 
the scissors operator approximation describes nicely the low lying excited 
states, which is seen in the good determination of the static dielectric func- 
tion and the low energy structures, i.e. E\ and E2, in the imaginary part of 
the dielectric function, it seems to fail for the description of the higher ex- 
cited states. This is not surprising because the higher excited states which are 
free electrons like are most probably well described within LDA and need no 
scissors-operator shift. This is supported by the fact that the the energy-loss 
function, -Ime -1 , within the LDA has it maximum roughly at the free electron 
plasmon frequency whereas within the scissors approximation its maximum is 
shifted to higher energies. For our purpose the scissors-operator shift remains 
a good approximation for the description of the low-lying excited states of 
semiconductors and their optical properties. 



6 Applications 

6.1 Optical Properties 

We have used our FP-LMTO method and the formalism outlined above to 
calculate the optical properties of materials H[lO 11,1^]. In general our results 



are often in good agreement with the experimental results. For semiconduc- 
tors, however, good agreement with experiment is only achieved when the so 
called scissors-operator shift is used. Figure 1 presents our relativistic calcula- 
tion of the imaginary part of the dielectric function of GaAs compared to the 
experimental results of Ref . J39| . The LDA relativistic results underestimates 
the band gap by about 1.3 eV. When the imaginary part of the dielectric 
function is shifted to higher energies by 1.3 eV the results the E\ and E2 
peaks are overestimated in our calculation. One needs to shift the spectrum 
by less than the band gap as done in Ref. jl(| to produce good agreement 
with experiment. It seems then that the optical band gap is less than the 
band energy gap (1.5 eV). The optical band gap is produced by interband 
transitions to the low lying conduction states. Excitonic effects are therefore 
important and are responsible for the reduction of the energy gap of semicon- 
ductors. It is interesting to notice though that the static dielectric function 
are in good agreement with experiment for GaAs, Si, and Ge when the shift 
correspond to the energy band gap obtained from photoluminescence PJTp|. 

More interesting are the wide band-gap materials where the LDA calcu- 
lated static dielectric function is in good agreement with experiment despite 
that the band gap is still underestimated by LDA. Correcting the band gap 
using the scissors operator makes the static dielectric much small than the 
measured value. As an example of wide gap material, we present in Figure 2 
and 3 the imaginary part of the dielectric function of GaN for the cubic (3) 
and wurtzite structure (B4). 

Table I shows that our LDA dielectric constant calculations are in agree- 
ment with available experimental results and the pseudo-potential (PP) re- 
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6j of GaAs 
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Fig. 1. Calculated Imaginary part of the dielectric function of GaAs at the experi- 
mental equilibrium volume both within LDA and shifted by 1.3 eV, compared with 
the experimental results of Ref. The experimental E\ is only slightly under- 
estimated while Ei is overestimated. Notice that the shifted dielectric function by 
1.3 eV, which produces the correct band gap, overestimates the peak positions by 
about 0.3 eV. Excitonic effect should shift these peaks to lower values in agreement 
with experiment. 

Table 1. Calculated static dielectric function taa for GaN compared to pseudopo- 
tential (PP) results and experiment. For the wurtzite structure we have calculated 
eL for a polarization parallel to the xy plan and e^ which is perpendicular. 





zinc-blende 

e x> 


wurtzite 

c oc 




PP 


5.74 


5.48 


5.60 


Present work 


5.96 


5.54 


5.65 


Expt. 




5.35 


5.35±0.2 



suits Q including local-field effects (an error about our calculation is re- 
ported in Ref. our value for eL is not 4.48 but 5.54 and the PP value 
should then be 4.48). It is interesting to notice that static dielectric is in good 
agreement for the for all the nitrides [[l0] while the band gap is underesti- 
mated. The scissors-operator shift fails to explain the static dielectric function 
of large gap semiconductor. Recently, both local-field effects and electron-hole 
interaction were included on an ab-initio computation of the dielectric func- 
tion of few semiconductors pl| , ^2| by extending the semi-empirical Hanke 
and coworkers approach p^ , [43| which is based on the solution of the Bethe- 
Salpeter equation G3. The excitonic effects seem to improve significantly 
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Fig. 2. Calculated imaginary and real parts of the dielectric function of GaN in its 
cubic and wurtzite forms. The LDA band gap of the cubic phase is 1.8 eV and the 
wurtzite phase is 2.2 eV. 



the agreement between theory and experiment. However for large band-gap 
semiconductors, such as diamond, the inclusion of the excitonic effects seem 
to underestimate the optical band gap by as much as 1 eV Ell . It is not clear 
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from these calculations whether the static dielectric function for wide-band 
gap semiconductors is improved when excitonic effects are included. More 
theoretical work along these lines is needed to fully understand the dielectric 
function of wide-gap semiconductors. 

6.2 Magnetic Circular Magnetic Dichroism 

X-ray absorption spectroscopy (XAS) probes selectively each core orbital of 
each atomic species in a material. Two decades ago the theoretical work of Er- 
skinc and Stern show that the x-ray absorption could be used to determine the 
x-ray magnetic circular dichroism (XMCD) in transition metals when left and 
right circularly polarized x-ray beams are used [H . More recently these ideas 
were implemented experimentally and XAS was used to determine the local 
magnetic properties of each magnetic atomic orbital in a magnetic compound 
p5jE(|. Thus the circular magnetic x-ray dichroism is an important tool for 
the investigation of magnetic materials gl^g^loljl^QIIIIJIIQ, 
especially through the use of sum rules for the direct determination of the 
local orbital and spin contributions to the total magnetic moment [ ^0[ . 

Thole and co-workers show that the circular-magnetic-x-ray dichroism is 
related to the magnetic moment of the photo-excited atom when the core elec- 
tron is excited to the conduction states that are responsible of the magnetic 
properties of the material. On the theoretical side, Ebert and his co-workers 
[ pUpi^ ] have developed a fully-relativistic local-spin-density-approximation 
approach that was used with success to calculate the XMCD at the K-edge 
of Fe, the L 3 -edge of Gadolinium, and Fe and Co multilayers. Wu et al used 
slab linear augmented plane wave method to study the 1^,3 XMCD of Fe 
[ j56"l . Brouder and co-workers uses Multiple-scattering theory to solve the 
Schrodinger using spherical potentials and spin-orbit coupling as a perturba- 
tion in the final state . Recently Ankudinov and Rehr used a method based 
on a non-relativistic treatment of propagation based on high order multiple 
scattering theory and spinor-relativistic Dirac-Fock treatment of the dipolc 
matrix elements to calculate the Fe K edge and Gd L 3 edge XMCD Q . 

The calculation of the x-ray absorption for left and right circularly po- 
larized x-ray beams is implemented within the local-density approximation 
(LDA) by means of all-electron full-relativistic and spin-polarized full-potential 
linear muffin-tin orbital method (LMTO). The core electrons are spin-polarized 
and their electronic states are obtained by solving the full-Dirac equation, 
whereas for the valence electrons the spin-orbit coupling is added pertur- 
batively to the the semi-relativistic Hamiltonian. The total Hamiltonian is 
then solved self-consistently. To calculate the polarization dependent cross- 
section we consider the case where the internal field polarizes the spins along 
the magnetization easy axis. With respect to this axis we defined the left- 
and right-circular polarization, which correspond to the photon helicity (+h) 
(-U) respectively and the following dipole interaction: e±p = -j=(\/ x ±iV y ). 
The absorption cross-section fi± for left (+) and right (— ) circular polarized 
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Fig. 3. Calculated x-ray absorption at the K-edge of Fe for left and right circularly 
polarized light compared to the experimental spectrum. The difference between the 
two spectra (barely visible on the graph) represents the x-ray magnetic circular 
dichroism. 



x-ray calculated at the relativistic j± (£ ± |) core level and in the dipole 
approximation is given by: 

M±M = ^ ^2(j±m J± \e ± p\nk)(nk\pe ± \j ± m j± )5(uj~E nk +E J± )(36) 

using LDA in conjunction with the relativistic full-potential LMTO tech- 
nique. 

Figure 3. represent the K-edge x-ray absorption of Fe, for left and right 
circularly polarized light, compared to the experimental results. The agree- 
ment at low energy with experiment is good and start degrading at higher 
energies above the mean absorption peak. It is of interest to point out that 
the magnetic x-ray dichroism at the K-edge which is due to the spin polariza- 
tion and the spin-orbit in the final state is very small in the case of Fe. The 
difference between the right and left circularly polarization of the light is not 
even visible on the graph. However, the x-ray magnetic circular dichroism can 
be measured and Figure 4 shows a good agreement of the calculated dichroic 



signal with the experimental results of Shiitz 1 46 . 

At the L/2,3 edge of 3d transition metals the x-ray magnetic dichroism 
is much important because it is meanly due to the presence of the strong 
spin-orbit coupling in the initial 2p states (in the case of Fe the spin-orbit 
splitting between the 2p 3 / 2 and 2p 1 / 2 is about 13 eV). In Figure 5 we show 
the calculated x-ray absorption and XMCD at the Co in PtCo ordered alloy 
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Fig. 4. Calculated x-ray magnetic x-ray dichroism at the K-edge f Fe compared to 
the experimental spectrum of ShutzJ4(J. 
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Fig. 5. Calculated x-ray absorption and magnetic x-ray dichroism at the L2,3-edge 
of Fe compared to the experimental spectrum of Grange et al. |2lj . 



To compare the results with experiment we have to take into account the 
effect of the core hole and the experimental resolution. This is done by convo- 
luting the calculated spectra by a Lorentzian of widths of 0.9 eV and 1.4 eV 
for the L2 and L3 edges, respectively, in addition a Gaussian broadening of 
0.4 eV is added to take into account the experimental resolution. The calcu- 
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Fig. 6. Calculated x-ray magnetic x-ray dichroism at the L2,3-edge of Pt compared 
to the experimental spectrum of Grange et oZ.0|. 

lation of the x-ray magnetic circular dichroic signal ignoring the electron-hole 
recombination effect provides a semi-quantitative agreement with the exper- 
imental spectra. Hence, we believe that the core hole effect represented here 
by a Lorentzian broadening plays a significant role in determining the correct 
L 3 /L 2 branching ratio for 3d transition metals. The underestimation of the 
L2.3 branching ratio remains a challenge for theorists and further theoretical 
development along the line proposed by Schwitalla and Ebert (5^] is needed 
to bring the theory at the level of the experiment. 

For the 4d-transition metals, the core hole is deeper, and the agreement 
with experiment of the XMCD is satisfactory. Figure 6 shows the calculated 
XMCD at the site of Pt of the CoPt ordered alloy. 

In contrast to what is obtained for Co, the results for the Pt site show a 
much better agreement with experiment, due to the fact that the core hole 
effect is less intense (core hole much deeper than that of Co) . For the Pt atom 
we used both a Lorentzian (1 eV) and a Gaussian (1 eV) to represent the core 
hole effect and a Gaussian of 1 eV width for the experimental resolution. The 
experimental and theoretical Li and L3 edges are separated by a spin-orbit 
splitting of the 2p core states of 1709 and 1727 eV respectively. The width 
of both L2 and L3 edges is comparable to experiment, but the calculated L2 
edge is much larger. This produces a calculated integrated branching ratio of 
1.49 which is much smaller than the experimental ratio of 2.66. Here again 
the theory is underestimating the branching ratio. 

7 Conclusion 

We have reviewed the FP-LMTO method and the implementation of the 
optical properties and x-ray magnetic dichroism within the local density ap- 
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proximation. We have showed that the momentum matrix elements can be 
evaluated as a muffin-tin contribution and a surface term. The method has 
been successfully used to compute the optical properties of metals ]lT| , p^ |, 
semiconductors p[|l(| and magneto-optical properties |Tl"| o f transition met- 
als alloys, as well as x-ray magnetic circular dichroism |2C| , pl| , p2[ with high 
precision. 

For small-gap semiconductors a scissors-operator shift should be used to 
reproduce the static and dynamic dielectric function JlO| ]. Excitonic effects 
seem to be important in reproducing the correct optical energy gap f4l"| , [42|| . 
For wide-gap semiconductors the local-density approximation (LDA) static 
dielectric function is in good agreement with experiment and no scissors- 
operator shift is required despite the underestimation of the band gap by 
LDA @. 

For the computation of the x-ray magnetic circular dichroism the agree- 
ment with experiment is rather good [)20[|2 l|f^ , [in] , p r 2T|^ , |^ , [55f . However, the 
so called branching ratio is underestimated by the theory. More theoretical 
work where the electron core-hole interaction is taken into account is needed 
to bring the theory at the quality level of experiment f37f . 

Part of this work was done while one of us (M.A) was at Ohio State 
University and were supported by NSF, grant number DMR-9520319. Su- 
percomputer time was granted by CNUSC (project geml917) on the IBM 
SP2 and by the Universite Louis Pasteur de Strasbourg on the SGI O2000 
supercomputer . 
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